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Abstract 

A continuum of monopole, dyon and black hole solutions exist in the Einstein- 
Yang- Mills theory in asymptotically anti-de Sitter space. Their structure is studied 
in detail. The solutions are classified by non-Abelian electric and magnetic charges 
and the ADM mass. The stability of the solutions which have no node in non- 
Abelian magnetic fields is established. There exist critical spacetime solutions which 
terminate at a finite radius, and have universal behavior. The moduli space of the 
solutions exhibits a fractal structure as the cosmological constant approaches zero. 
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1 Introduction 

For a long time, it was believed that no regular particle-like stable solutions (solitons) with 
finite mass can exist in self gravitating systems unless the stability is guaranteed topo- 
logically. The Einstein theory in vacuum and the Einstein-Maxwell system do not admit 
solitons. It came as quite a surprise when Bartnik and McKinnon (BK) found globally 
regular solutions to the SU(2) Einstein Yang-Mills (EYM) theory without scalar fields 
It was unexpected to find that self gravitating Yang-Mills systems produced solitons. Un- 
fortunately, the BK solutions were shown to be unstable against linear perturbations 0. 
Later, other fields such as Higgs scalar fields and dilaton fields were included in the EYM 
action, but with the exception of the Skyrmions, all turned out to be unstable (see || for 



Interest in the BK solutions was renewed with the discovery of black hole solutions to 
the EYM equations @, §. These non-Abelian black holes aparently violate the no-hair 
conjecture ||. But these non-Abelian black hole solutions are also unstable, and again 
other fields were added in the hope of achieving stability without success (see Ref.'s 0, 0]) 
for a review). 

We stress that it is a surprise that there are static solutions to the Einstein Yang- 
Mills equations at all. There are no static solutions to the Yang-Mills equations in 4 
dimensional flat space. We can see this with a simple argument given by Deser ||. The 
conservation of the canonical energy momentum tensor, d u T^ v = 0, implies that for a static 
field configuration djT^ = 0. The total divergence of the quantity x l Ti 3 must vanish to 



maintain finite energy and regularity, / d d l xdj{x l T i 3 ) = 0. But dj(x l Ti J ) = Tj J +x l djT i J = 



Since the integrand above is positive definite for d = 4, Fij and F 0i must vanish. Thus 
there are no regular static solutions. 



The argument above cannot be extended to curved spacetime. The conservation law 
T^y, = leads to 



a review 



Tj J so that 




(1) 
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The failure of Deser's simple argument in curved space implies the possibility of having 
static solutions in curved space. Gravity supplies the attractive force needed to balance 
the repulsive force of Yang-Mills gauge interactions. Indeed, any solution to SU(2) EYM 
equations in asymptotically Minkowski space which is regular asymptotically is also regular 
for all r > g]. 

The particle-like and black hole solutions were later studied in a cosmo logical context. 
The behavior of static solutions to the Einstein Yang-Mills equations depends considerably 
on the sign of the cosmological constant. The solutions can be separated into two families; 
A > and A < 0. The solutions where A = are the BK solutions. Their asymptotically 
de Sitter analogs (A > 0) were discovered independently by Volkov et. al. and Torii et. 
al. fllU|. The BK solutions and the cosmological extensions all share similar behavior, 



and are unstable |TT], |12|| . (See Ref. || for a review). Recently, asymptotically anti-de 
Sitter black hole solutions [|l3j and soliton solutions fl4|, [15| were found which are strikingly 
different from the BK type solutions. In particular, the asymptotically anti-de Sitter AdS 
EYM equations have solutions where the field strengths are non-zero everywhere. These 
solutions were also shown to be stable against spherically symmetric linear perturbations. 
These solutions are the only EYM solutions solutions that are stable. This discovery would 
be very important to cosmology if the universe was ever in a phase where the cosmological 
constant is negative. 

Another new feature of the EYM theory in AdS is the existence of dyon solutions. If 



A > the electric part of the gauge fields is forbidden |16] if the ADM mass is to remain 
finite. Scalar fields must be added to the theory in order for the boundary conditions at 
infinity to permit the electric fields and maintain a finite ADM mass JI5|, [0|. 

Recently a tremendous amount of interest has evolved in field theories in AdS space. 
There is the AdS/CFT correspondence [17|. Conformal field theories in d dimensions {Rd) 
are described in terms of supergravity or string theory on the product space of AdS^+i and 
a compact manifold. There are intimate relations between data on the boundary Rd of 
AdS^+i and data in the bulk AdS^+i- In the present paper we are examining the Einstein- 
Yang-Mills theory in asymptotically AdS space. The boundary in space must be playing a 
crucial role for the existence of stable monopole and dyon solutions, more detailed analysis 
of which is, however, left for future investigation. We also note that in the three-dimensional 
AdS space there exist nontrivial black holes p0[ and monopole/instanton solutions [ f2l|j . 
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When the value of the cosmological constant A is varied, the space of monopole and 
dyon solutions, the moduli space, also changes. With a finite negative A, solutions exist in 
continuum. They are classified in a finite number of families, or branches. With a vanishing 
or positive A solution exists only in a discrete set, but there are infinitely many. One natural 
question emerging is how these finite number of branches of solutions in continuum become 
infinitely many discrete points as A < approaches 0. There is a surprising hidden feature 
in this limit. We shall find a fractal structure in the moduli space, which seems to explain 
the transition. 

In the next section the general formalism is given and the equations of motion are 
derived with a spherically symmetric ansatz. Conserved charges in the Yang-Mills theory 
is defined in section 3. Some general no-go theorems are derived from sum rules in Section 4. 
New soliton solution in asymptotically anti-de Sitter space are explained in section 5. The 
critical spacetime which have universality near the edge of the space is also examined. Black 
hole solutions which have both magnetic and electric non-Abelian charges are presented 
in Section 6. The dependence of the moduli space on the cosmological constant A is 
investigated in Section 7 where the fractal structure is revealed when A approaches zero 
from the negative side. The detailed analysis of the stability of the monopole solutions is 
presented in Section 8. The subtle boundary condition in the problem requires elaboration 
of the previous argument presented in the A = and A > cases. 



2 General Formalism 

In non-Abelian gauge theory, the field equations have solutions which exhibit a magnetic 
charge. In the 't Hooft-Polyakov monopole solution 



<&a = — H(er), A° a = 0, A\ = -e aij ^(l - K(er)). (3) 
er er 



where $ a is a triplet Higgs scalar field. Its stability is guaranteed by the topology of the 



triplet Higgs scalar field. [22] The U(l) magnetic charge takes a quantized value, Air/e. 
Dyon solutions were obtained |23| starting with the above ansatz (|3]) but with a non-zero 
value for A®, (i.e. A° a = (x a /er 2 )J(er)). 
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In this paper we look for monopole and dyon solutions in the Einstein- Yang-Mills theory 
without scalar fields; 



S = J d A x^f^g 



—^—(R - 2 A) - -F atxv F a uv 



(4) 



The Einstein and Yang-Mills equations are given by 

W v - \g» v {R - 2 A) = 8nG 

We suspect that the gravity provides attractive force to balance the equation. 
We look for spherically symmetric solutions. The metric takes the form 



(5) 



ds 2 



H dr 2 

— dt 2 + + r 2 {d6 2 + sin 2 9 dcj) 2 ) 



P 



H 



(6) 



whereas Yang-Mills fields are given [24, |25], in the regular gauge, by 



A (o) 



^-{A ^dt + A 1 ^dx k 
2e r r z 



01 / XjX k \ 1 - <> 
H Ojk T~ ) ax k — Cjkl T 

1" \ T** J T* 



-x k dxi 



(7) 



Here the Cartesian coordinate x fc 's are related to the polar coordinates (r, 0, 0) as in the 
flat space. H, p, A Q , A x , 0i and 02 are functions of r for monopole or dyon solutions. In 
the discussion of the stability of the solutions they depend on both t and r. The regularity 
of solutions at the origin demands that H, p are finite, whereas Aq, A\, and 0i — > and 
2 — * 1 at r = 0. 

2.1 Simplification of the static gauge field ansatz 

Let A = A^dx^ = |rM."<i;r M , where r a are the usual Pauli matrices. In terms of the 
basis in spherical coordinates (r r ,rg,T^) 
(i = r,9,(f)), the ansatz (]?]) is written as 

1 



(n r , no, n^)r which satisfies [t^t,] = He^Tk 



A {0) = — \A T r dt + A X T r dr + (0it + (0 2 - l)r (f> )d9 



+ (-(02 - 1)t0 + 0ir ) sin 



Note that there are no singularities in this gauge. Next make a gauge transformation 
A = SA^S- 1 - (i/e)dS ■ S- 1 where 



/ _|_ e i(0+C)/2 ^ . ' <>••>• V \ 



s 



cos - + e-^'^ 12 sin - 
2 2 



V 



_ e W-n)/2 sin B _ + e -i(<p+n)/2 cog °_ 
2 2 / 



, n = n(t,r) . 



Useful identities are 



Sr r S- 1 
StqS -1 

2idS ■ S- 1 



cos VLt\ — sin Qt2 
sin Qti + cos Qt2 

— (fl'dr + J7c?t)r3 — G?0(sinfb~i + cosf^) 
+d0(sin 6 1 cos f)ri — sin sin Qr 2 — cos #r 3 ) 



(9) 



(10) 



The new gauge potential is 

A = — VuTzdt + vT^dr + (wt\ + wr 2 )d9 + (cot #r 3 + wt 2 — wti) sm9d(p\. (11) 

where 

it = A + Q 

v = Ax + Q,' 

w = +0i cos Q + 4> 2 sin f2 

w = — 0i sin Q + 02 cos f2 . (12) 

Note that the gauge transformation (g) is singular at 6 = and tt. Eq. (|T3) is the 
gauge potential in the singular gauge. It has a Dirac string. One can always choose 
Q(t, r = 0) = 7r/2 with which the boundary conditions at r = are u = v = w = and 
w — 1. With appropriate f2(t, r) one can set z/(t, r) = or u(t, r) = 0. 

A straightforward calculation leads to the field strength F = dA — ieA A A: 



— \ (z> — u')T 3 dt A dr + 
2e 



flO — UWjTi + [W + UW)T2 



dtAd6 



[UW + W)T\ + TO — W)T2 



dt A sin 



+ («/ — vw)ti + (u>' + wv)t2 dr A d# 
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+ (w' — vw)t2 + (—w — vw)t\ dr A sin 
-(1- w 2 -w 2 )r 3 d9 A sin 0d<j>\ . (13) 
The configurations where u — u — 0, w = w =constant, and w 2 + w 2 = 1 are pure gauge. 

2.2 Equations of motion 

In the general spherically symmetric metric @ tetrads are 



H 1 

eo = - — dt , ei = —j=dr , e 2 = rdfl , e 3 = rsin^d^ . (14) 
P V# 

In the tetrad basis F a t, = {e a ) IJi (eb) 1/ F' lu and the energy-momentum tensors are = 
F^F^ — F de W. The nonvanishing components of the Yang-Mills equations ([5]) 



arc 



(pr 2 (u — z>)) — 2-^{w(to + ty) + w(uw — w)} = 
H 



(2 / / \ \ f / / / 2 ~2\~i 

pr (u — u)j ^ — 2 — { — + w u; + + w )} = 



— (w' — wu)\ — ( — (w — uw 
V J \H 



t 



P , ~n w(l-w 2 -w 2 ) H , , 

+— « ™ + to J H v[w + wz/j = (J 

H pr 2 p 

— (w' + wu)) — ( —(w + uw* 
v J \H /./ 



V , ~ • \ w(l — w 2 — w ) H 

+—u(uw — w)-\ 1 v(w — wv) = (J . (15) 

H pr 2 p 



The nonvanishing components of the energy-momentum tensor are given by 

Too = \(A + B) 
T u = - 2 {-A + B) 

T 2 2 = T 33 = —A 

e 2 
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T i = -\c (16) 



where 



1 1 

A 2/ ■ /\ 2 i / 1 2 ~2\2 



v 2 H 
B = -^—{(uw + w) 2 + (w - uw) 2 } + — {(w' - vwf + (w' + vw) 2 } 



C = + vw){uw + w) + (w f — vw){w — uw)} . (17) 



r 2 H 
2p 

r 

The Einstein equations reduce to 

p' 8ttG rB 
_ _ __ 



p e 2 H 

H' l-H 8ttG , . . 
+ 5— = — =- (A + 5) + A 

(18) 

It is convenient to introduce m(r) defined by 

, , 2m(r) Ar 2 , . 

H(r) = l —-^r- • 19 

r 6 

m(r) is the mass contained inside the radius r. p(r) =constant and m(r) = corresponds 
to the Minkowski, de Sitter, or anti-de Sitter space for A = 0, > 0, or < 0, respectively. 
Then the second equation in ( pf) becomes 

m ' = ^ r'HA + B) . (20) 

The system of the Einstein- Yang-Mills equations contains one redundant equation. The 
third equation in ( ]T8|) follows from ([T5|) and the rest of fll8]) . 

2.3 Static configurations 

It is most convenient to take the v — gauge for static configurations. The second 
equation in ( |15D then yields ww' — w'w = 0, which leads to w(r) = Cw{r). By a further 
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global rotation Q =constant in ( |T2] ) one can set w = 0. As a result 
A = — ^UT^dt + wT\dQ + (cot 6t 3 + wr 2 ) sin ( 



F = — | — u'r 3 dt Adr + uwdt A (r 2 <i# — T\ sin 



+w'dr A (rxd9 + t 2 sin 9d(f)) - (1 - w 2 )r 3 d8 A sin 
Then the Einstein- Yang-Mills equations are 



/ 2 /V ^P 

(r pu) = — 



P 2 

-uw 
2p 2 



2i- 



P 



m 



-P 



u> (1 — w 2 ) 



, 2 U 2 W 2 p 2 



H 2 



^^+\r 2 Aur + H( W r + u ^ 



(21) 

(22) 
(23) 
(24) 

(25) 



where v = 4irG/e 2 . 

These equations are solved with the given boundary conditions. Near the origin solu- 
tions must be regular so that 



u(r) 
w(r) 
m(r) 
p(r) 



ar + -{ - 2b + -A + 2v(a 2 + 46 2 )}r 3 

5 3 
1 - br 2 



-via- + 46 )r 



2\„3 



v(a 2 + Ab 2 )r 2 



(26) 



where a and b are arbitrary constants. The boundary conditions at the origin of the EYM 
equations are completely determined by the values of the constants a and b. 

At space infinity the energy-momentum tensors T a ^ in ( pp must approach zero suffi- 
ciently fast. Further we expect that the metric must asymptotically (anti-) de Sitter space, 
depending on the value of A. This, with the equations of motion, leads to the asymptotic 
expansion at large r; 



u = u + Mi- H 

r 

m = M + mi — I - ' ' 

r 



w = w + W\~ + 
r 

1 

P = P0+P4^ + 



(27) 



where u , u\, Wq, Wi, mi, po and p 4 are constants to be determined and M is the ADM 
mass, M = m(oo) — m(0). 
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3 Conserved charges 



Solutions to eq.'s ( p2|) to (p5j) are classified by the ADM mass, M = m(oo) — m(0), electric 
and magnetic charges, Qe and Qm- From the Gauss flux theorem 



' e \ e 
.Qm) ~ 4^ 



(28) 



are conserved. With the ansatz in the singular gauge (|TT|) and the asympotitic behavior 
(p7|)) the charges are given by 

w»J u-»„ 2 y 2 

Notice that the electric charge Qe is determined by u±, whereas the magnetic charge Qm 
by Wq. If (u,w,m,p) is a solution, then (—u,w,m,p) is also a solution. Dyon solutions 
come in a pair with (±Qe, Qm, M). 

The charges are not gauge invariant, however. Under a local gauge transformation 
v4 — > UAU^ 1 — (i/e)dUU~ 1 , Qe and Qm are transformed to 



jdS k ^~gU{x) ( \u-\x) (30) 



In non-Abelian gauge theory a set of charges {Qe U , Qm U } are conserved. In the rest of 



the paper we use the charges, , defined in the singular gauge. 

The effective charge Q e r |l] is defined by the asymptotic behavior of H(r); 

H(r) = l- — + %-^Ar 2 . (31) 
r r z 3 

In terms of the coefficients in (|27|), Q^ s = —2m\. This requires that m\ < which indeed 
is the case. After inserting eq. (|27|) into eq.'s (|22|) and we find the relation 



Q e 2 fr = 2t;Tr(Q| + Q| f )-^^i. (32) 

3 Po 



Eq. fl24j ) implies that p(r) is a monotonically decreasing positive function so that po > 
and p4 > 0. The effective charge is smaller (larger) than 2vTr(Q E + Q 2 M ) for A > (< 0). 
The relation fl32f) incidentally implies that the charges defined in the singular gauge have 
physical, gauge invariant meaning. 
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4 Sum rules 



Sum rules are obtained from the equations of motion. First, multiply both sides of fl23|) by 
u and integrate in part. 



pr 2 uv! 



= £ 2 dr^r 2 p(u') 2 + 2j^u 2 w 2 ^j . 



(33) 



Secondly, multiply both sides of ( p2|) by w and integrate in part: 



H , 

— WW 



P 



f d r l-(w') 2 -^V-Vfl-w 2 ) 1 • (34) 
Jn { p H pr z j 



n 

Thirdly, divide both sides of Q22I) by w and integrate in part: 



Hw' 



pw 



ri 



p 2 H ( w'\ 2 1 



dr + - ~ • (35) 



ri 



p \w J pr 



These relations are valid, provided the integrals on the right hand sides are defined. Several 
important conclusions follow from (^) - (|35|) . 



4.1 In asymptotically flat space 

Consider fl33|) with ri = and r 2 = oo. For regular solutions u(0) = 0. Both p and H 
approach constant as r — ■> oo. The finiteness of the ADM mass requires that uw\ r=OQ = 0. 
In the expansion u w = 0. On the other hand, if u ^ 0, w = and Eq. (p3|) implies 
lOi 7^ so that Eq. (p^ ) cannot be satisfied. Hence m(oo) = 0. Then the left hand side of 
(|33|) vanishes, implying that u(r) must vanish identically. There is no regular electrically 
charged solution. Furthermore, (B^) can be solved only if (w ) 2 = 1 as H(oo) = 1, therefore 
the magnetic charge Qm vanishes. 

Suppose that w(r) never vanishes and w 2 < 1 for < r < oo. Consider fl3"5]) with 



r\ = and r 2 = oo. The l.h.s. vanishes, but the integrand on the r.h.s. is positive definite 
except for the pure gauge configuration w(r) = ±1. This implies that non-trivial solutions 
with w 2 < 1 must vanish at least once. 

We also note that the singular solution w(r) = 0, u'(r) = r~ 2 , and p(r) = 1 is nothing 
but the Reissner-Nordstrom solution. 
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4.2 In asymptotically de Sitter space 

In asymptotically de Sitter space H(r) — > — Ar 2 /3 as r — > oo. In this case the finiteness 
of the ADM mass does not forbid non- vanishing uw at r = oo. However, there arises a 
cosmo logical horizon at r = where H(rh) = 0. 

It follows from (^3) and fl25|) that u or w must vanish at r = r^. Now consider 
(|35|). Suppose that w(rh) ^ 0, or equivalently w ^ for ri < r < r 2 = for some 
ri > 0, the left hand side is finite so that u{rh) = 0. However, Eq. fl23|) implies that 
u(r) is a monotonically increasing or decreasing function. With the boundary condition 
u(0) = 0, the only possibility available is u(r) = 0. This conclusion remains valid even if 
w(rh) = 0. In this case the first and second terms on the r.h.s. give positively divergent 
contributions near r^, whereas the l.h.s. remains finite. To summarize, there is no solution 
with nonvanishing u(r). 

If w > for < r < rh, then the left hand side vanishes in the r\ — > and r 2 — > 
limit. If one further assumes that w 2 < 1 in the interval, the integrand on the right hand 
side is positive definite so that the only solution is w(r) = ±1, which is a pure gauge. 

This argument also shows that a nontrivial solution w(r), which satisfies w 2 < 1 for 
< r < r h , must vanish at least once in this interval. We have numerically looked for 
solutions in which w 2 > 1 for small r, but have found that no such solution exists. 

4.3 In asymptotically anti-de Sitter space 

In asymptotically anti-de Sitter space there exist solutions in which H(r) > everywhere. 
As H(r) ~ |A|r 2 /3 for large r, the condition for the finiteness of the ADM mass requires 
only that uw — > constant. In other words, both u and w may approach nonvanishing 
values as r — > oo. Consequently, with the expansion (^) the l.h.s. of ( p3|) with r x = and 
r 2 = oo is — PqUqUi and can be non-vanishing. Solutions with u(r) ^ are allowed. 

In critical cosmological horizon appears and H(rh) = 0. In this case the argu- 

ment above for the asymptotically de Sitter case applies and either u or w must vanish at 
r = We shall find, indeed, wijh) = below. 

5 Soliton solutions 
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5.1 The BK solution 



Particle like solutions of the EYM equations in asymptotically Minkowski space were first 
found by Bartnik and McKinnon (BK) []J in 1988. When u = 0, we can solve the equations 
( p2|) to (^) numerically with the boundary condition a = 0. 

As already discussed, w(r) = 1 and H(r) = 1 corresponds to a pure gauge configuration. 
Similarly, if w = as r — > oo we are left with the RN solution which is singular at r = 0. 
The Yang-Mills charge (E^) vanishes in asymptotically Minkowski space since wq = ±1. 
The effective charge Q c s = in all BK solutions, since Qm = A = 0. All solutions possess 
a metric that is asymptotically Schwarzschild. 

There is a discrete set of BK solutions, labeled by the number of nodes in w, n 6 [1, oo), 
and the free shooting parameter b. Since all BK solutions have at least one node, n > 1, 
the solutions are unstable against spherically symmetric perturbations [§]. 



5.2 Solutions with a cosmological horizon 

Solutions in asymptotically de Sitter space were obtained by adding a cosmological constant 
(A) term to the Einstein equations. Solutions display the same basic properties as the BK 



solutions |T0[. Just as in the BK solution, these equations are solved numerically, using the 
shooting method and requiring a = 0. 

In asymptotically de Sitter space, a cosmological horizon, where H = 0, develops at 
r = Th- At the horizon H'(rh) ^ 0. Near the horizon 

w{r) = w + WiX + w 2 x 2 + ■ ■ ■ 

p(r) = p + Pix + p 2 x 2 H 

m(r) = m + m^x + m 2 x 2 + ■ ■ ■ 

H(r) = h 1 x + h 2 x 2 + --- (36) 

where x = r — r h . m and r h are related by 1 — (2m /r h ) — (Ar^/3) = 0. With given w , 
r h , and p , the equations (|22|), (|24D , and (|^) (with u = 0) determine all other coefficients, 
provided w ^ 0, ±1 and p ^ 0; 



v 



.2 i \2 



h 
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2m\ 2m 2 

hi = 1 n o Ar h 

r h ri 3 



Wi 



w (l - wl 
h\r\ 



Pi = 37 

r h 

and so on. This expansion is valid independent of the value of A. The critical case 
hi = H'{rh) = requires a special treatment, and will be analyzed in Section (5.3). 



Kiinzle and Masood-ul-Alam || have argued that w(r) has y\r — singularity at r = r^. 
However, we have found that the regular expansion (^) is valid. 

Just like the BK case, there are a discrete set of solutions labeled by the number of 
nodes in w(r), n and the parameter b. Solutions in w and m, have the same form as the 
BK solutions except that w(r) no longer approaches 1 at infinity. Eq. ([32|) implies that 
there is a non-vanishing charge Qm (or w ^ 1) if A ^ and that the solutions are all 
asymptotically Reissner- Nordstrom type. |w(oo)| is slightly greater than 1 for the n = 1 
solutions. When n > 2, 1 > w(oo) > where w(oo) — > as n — > oo. 

The mass m(r) also stays finite. For all indices n, it is small near the origin and does 
not grow until it approaches the horizon where it quickly climbs to a value near 1. After 
the horizon it stays almost constant. 

The position of the horizon depends on A. As long A is below some critical value, 
the geometry approaches Reissner-Nordstrom-de Sitter space in the asymptotic region as 
indicated by Eq. (0). Above some critical value the topology changes and a singularity 
appears. Since the position of the horizon depends on A, there is a value for A where 
this singularity and the horizon meet, in which the topology becomes a completely regular 
manifold. Above this value for A, the solutions are no longer regular. More details of the 



topology dependence on A can be found in JlO|- All the solutions are unstable [0. 



5.3 Solutions in asymptotically anti-de Sitter space 

As already discussed, there are no boundary conditions that forbid a solution to the EYM 
equations in asymptotically anti-de Sitter space with a non-zero electric component, u(r), 
to the Yang-Mills fields. Solutions to Eqs. (22) to (|25|) are determined with the cosmological 
constant A fixed at some negative value. 
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(a) Monopole solutions 

Monopole solutions are obtained by setting a — (u — 0). By varying the initial 
condition parameter b, a continuum of monopole solutions are found which are regular in 
the entire space. Just as in the BK and dS solutions, w crosses the axis an arbitrary number 
of times depending on the value of the adjustable shooting parameter b. In contrast to 
the A = and A > cases which have a discrete set of solutions in b and n, there is a 
continuum of solutions in b for each n. Typical solutions are displayed in fig. 1. 




Figure 1: Monopole solutions for A = -0.01 and v = 1. (a, b) = (0, 0.001) and (0, -0.001). 
(w,m) at r = oo are (0.318,0.035) and (2.304,0.201). 



The behavior of m and p is similar to that of the asymptotically dS solutions ||10|| . In 
contrast, as shown in fig. [I], there exist solutions where w has no nodes. These solutions 
are of particular interest because they are shown to be stable against linear perturbations, 
(b) Dyon solutions 

Dyon solutions to the EYM equations are determined if the adjustable shooting param- 
eter a is chosen to be non-zero for a given negative A. Just as in the monopole solutions, 
we find a continuum of solutions where w crosses the axis an arbitrary number of times 
depending on a and b. Also similar to the monopole solutions is the existence of solutions 
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where w does not cross the axis. As shown in Fig. ^ the electric component, u, of the EYM 
equations starts at zero and monotonically increases to some finite value. The behavior of 
w, m , H, and p is similar to that in the monopole solutions. 

Just as for the monopole case, dyon solutions are found for a continuous set of param- 
eters, a and b. For some values of a and b, solutions blow up, or the function H(r) crosses 
the axis and becomes negative. 




Figure 2: Dyon solutions for A = —0.01 and v — 1. In the top figure w(r) has no node 
(n = 0), whereas in the bottom figure it has one node (n = 1). 

(c) Critical solutions 

As the parameter b is increased, the minimum of H(r) hits zero from above, i.e. H(rh) = 
H'{rh) = 0. This constitutes a special case and needs careful examination. Numerical 
studies indicate that this happens in a finite range of the parameter a. The critical solution 
exists for both u(r) = and u(r) ^ cases. One example of solutions near the critical 
value [(a, b) = (0.01, 0.69)] is displayed in fig. 3. H{r) becomes very close to zero at r ~ 1. 
It has (Q E ,Q M ,M) ~ (0.015,0.998,0.995). 

When b = b c , w and p vanish at r = as well. As p(r^) = 0, p(r) = for r > 77,. The 
space ends at r = r^. There is universality at the critical point. 
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Figure 3: Dyon solution for A = -0.01, v = 1, a = 0.01 and b = 0.69. At a = 0.01, the 
critical value is b c = 0.7104. H almost hits the axis around r = 1. At b c , p(r) = for 
r > Th- The space ends at rv 

The numerical integration of the differential equations indicates that m(r) and u(r) are 
regular at r = r/j. The appropriate ansatz for the critical solutions with H(rh) = H'(rh) = 
is, for y = r h — r > 0, 



u 
w 

P 
H 

m 



u + uiy + u 2 y H 

y a [w + wiy + w 2 y 2 H } 

y p {po +Piy + P2y 2 h — } 

y^{ho + h iy + h 2 y 2 + ■■■} 
m + m x y + m 2 y 2 H 



(3f 



Eq. (^) implies that 7 = 2. When u{r) ^ 0, Eq. (^) leads to a = \. When u{r) = 0, Eq. 
([24]), instead, implies that a = |. Other relations obtained from eqs. (p2|) - ( 25 ) are 



Ah 



prim 



2wqUo 







mi 



VWo 

v 



(39) 



17 



The value of the index (3 is unconstrained when u = 0. However, if u 7^ 0, the consistency 
of eq. (|22l) , for instance, demands that 2(3 be an integer. The smallest value for (3 which 



satisfies the first relation in (|39|) is (3 = 2, as h > 0. We have confirmed this by numerical 
studies. The relation (|19|) further implies that 

2m A 2 
= 1 ri 

r h 3 " 

m x = -(Ar^-l). (40) 



From the two relations for mi, one in (|39|) and the other in (fMj), is determined as a 
function of v and A < 0: 



ll = 2jAf(Vi + 4«|A|-i) • mi; 



To summarize, the indices in (|38|) are given by 



2 

The coefficients mo, mi, and 



a = \ , /? = 7 = 2 . (42) 



w = , h = -3 (43) 

are all determined by u and A only. We are observing the universality in the behavior of 
the critical solutions. The coefficients uo and po depend on a or b as well. 
For small u|A| 1 

1 r- 1 



2 4 

777 , h o (44) 

They are all determined by v only. This universal behavior is clearly seen in the solution 
in fig. 3 which is very close to the critical one. 

The meaning of the critical spacetime is yet to be clarified. The space ends at r = r/,. 
It defines a spacetime with a boundary. 

(d) Spectrum of monopole and dyon solutions 

Monopole and dyon solutions permit non- vanishing charges Qm arid Qe, although there 
are solutions where Qm = and Qe 7^ or where Qe = but Qm 7^ 0. Non-zero charges 
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Figure 4: Mass M is plotted as a function of magnetic charge Qm for monopole solutions 
at A = —0.01 and v — 1. The number of nodes, n, in w(r) is also marked. In the lower 
branch Q M = 1 at b = 0.00168. 

Qm or Qe ensures that Q e s ^ (see Eq. (^) ) so that solutions are asymptotically of the 
AdS Reissner-Nordstrom type. 

In fig. H the mass M is plotted as a function of Qm for monopole solutions at A = —0.01 
and v — 1. The behavior of the solutions near tb = b c = 0.7104 needs more careful analysis. 

Dyon solutions are found in a good portion of the Qe-Qm plane. There are solutions 
with Q M = but Q E ^ 0. Although Q M = 0, i.e. to(oo) = ±1, w(r) ^ 0. In the shooting 
parameter space {a,b), these solutions correspond not exactly, but almost to a universal 
value for b ~ 0.0054. See fig. |[ More surprising is the fact that Qm takes a quantized 
value — (47r) _1 / 2 at b = 0.0061 independent of the value of a within numerical errors. We 
have not understood why it should be so. 

Solutions with no node in w(r) have special importance, as they are stable against 
small fluctuations. (See section 8.) In fig. |6] the spectrum of nodeless dyon solutions are 
presented in the parameter space (a, b). Notice that a must be small enough (a < 0.005) 
even for b < 0. 

(e) Dependence of the coupling v on the solutions 

The ADM mass M depends on the value of the coupling v = AirG/e 2 . As shown in 
fig. 0, M increases as v gets larger and decreases when v gets smaller. With fixed (a, b), 
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Figure 5: Magnetic charge Qm of dyon solutions as a function of the parameter b is plotted 
with various values of the parameter a. At b ~ 0.0061, Qm is independent of a, taking the 
quantized value — (47r) -1 / 2 = —0.2821 within numerical errors. 

roughly M oc v . The Yang-Mills fields u and w are roughly independent of v. 



6 Black hole solutions 

Not long after the BK solutions were discovered, black hole solutions were also found to 
be contained in the EYM equations 0] if different boundary conditions were used. These 
solutions generated a large amount of further study, as they apparently violate the no- 
hair conjecture H. Later, EYM black holes were studied in a cosmological context by 



including a positive cosmological constant pi. These black hole solutions share most of 



the properties as the soliton solutions, including their instability. 

Recently, purely magnetic black hole solutions were found in asymptotically anti-de Sit- 
ter space [JDJ . These solutions are drastically different from their asymptotically Minkowski 
or de Sitter counterparts. There are a continuum of solutions in terms of the adjustable 
shooting parameter that specifies the initial conditions at the horizon. Furthermore, there 
exist solutions that have no node in w and are in turn stable against spherically sym- 
metric linear perturbations. Here, we discuss the solutions found by Winstanley |l3j and 



also present new dyon black hole solutions. We also discuss the apparent shrinking of the 
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Figure 6: Spectrum of nodeless dyons. 

moduli space when the magnitude of A is decreased. Similar to the particle-like solutions 
already discussed, the moduli space becomes discrete in the A — > limit. 

6.1 Boundary conditions at the horizon 

Black hole solutions are obtained numerically by specifying the boundary conditions at the 
horizon and shooting for regular solutions w, u, m and p for < r < oo. The location 
of the horizon, r^, and the value p(r^) > can be arbitrarily chosen by scaling of t and r. 
We look for solutions in which H(r) > for r > r^. As H(rh) = but p(r^) ^ 0, Eqs. (p2|) 
- ( |25j ) require that either u or w vanishes at the horizon. A stronger condition is obtained 
from the sum rule fl35|) with r\ = and r 2 = oo. Its l.h.s. is finite so that «(r/J = on its 
r.h.s. Hence we are led to the expansion 

W = W + W\X + ■ ■ ■ 

U = UiX + ■ ■ • 

p = 1 + P\X + • • • 
H = h\X H 

m = mo + m\x + ■ • • (45) 
where x = r — r^. We have chosen p(rh) = 1 without loss of generality. 



21 



Mass Spectrum, A — —0.01 




Figure 7: A typical solution showing the dependence of v on the mass m(r). A = —0.01, 
a = 0.003 and b = 0.001. w(r) and u(r) do not have much dependence on v. 

There are two adjustable shooting parameters, (a, b) = (ui,wo). After inserting the 
ansatz into eq.'s fll9D to fl25|) we find 



r h Ar 3 h 

m Q — 

2 6 

h x = —(l-Arl-2m 1 ) 

•h 



Wi 



Pi 



w (l - wl) 



r h h l 










wlu\ 
hi 


I < 4 


1 

r h 





The asymptotic expansion at large r is the same as in (^). The ADM mass is given by 
M = m(oo). 
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6.2 New electrically and magnetically charged black hole solu- 
tions 



Just as the soliton solutions, purely magnetically charged black hole solutions are obtained 
by setting the adjustable parameter a to zero. The behavior of the solutions are similar 
to that of the solitons (see Ref. |13[ for more information). The number of nodes n in w 
can be 0, 1, 2, • ■ -. The black hole monopole spectrum of mass versus charge is displayed 
in fig. 0. It shows the spectrum for the n = and n = 1 arms. 
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Figure 8: Black hole mass vs. magnetic charge Qm for the n = and n = 1 arms and for 
different values of A. 



Solutions with both magnetic and electric charge are obtained by giving a a finite value. 
Dyon black hole solutions are similar to the monopole solutions except that u is nonzero. 
At the horizon u starts at zero and monotonically increases asymptotically to a finite value. 
H starts at one and quickly diverges, p starts at one and remains almost constant. Typical 
black hole dyon solutions are shown in Fig. [|. 

Again black hole dyon solutions with no node in w(r) are stable against small spher- 
ically symmetric perturbations. The spectrum of those nodeless black hole dyons in the 
parameter space (a, b) is plotted in fig. [II]. Notice the similarity between fig. ^| and fig. [ID. 
The nodeless solutions exist only for small a = u\ < 0.0055. b = wq must be around 1. 
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Black Hole Dyon Solutions, A — 



-0.01, v=1, r,= 1 




Figure 9: Typical dyon black hole solutions with no node and one node. H(r) and p(r) were 
not plotted for the sake of clarity. The top figure corresponds to (a, b) = (0.004,0.99998), 
whereas the bottom to (0.01,0.9998). 

7 Dependence on A — fractal structure 

The soliton and black hole solutions depend non-trivially on the value of the cosmological 
constant A. It has not been well understood why the continuum of solutions for negative 
A become a discrete set of solutions in the A — > limit, and remain discrete for all A > 0. 



Just as fig. H shows for the black hole solutions, fig. 4 and fig. |TT] shows the spectrum in 



mass vs. magnetic charge Qm plane for a give A. The width of each branch for a given A 



gets smaller as A approaches zero. Fig. [TT| indicates that as A — > 0, the branches collapse 
to one point, the BK solution, as the continuum of solutions vanishes. It is still unknown 
mathematically why and how this occurs. 

We would like to point out that there is a fractal structure in the moduli space of the 
solutions. This is most clearly seen in the parameter b v.s. mass M plot as displayed in fig. 
|l"2"|. As A becomes smaller, a new branch appears. The shape of branches has approximate 



self-similarity. Similarly, in fig. [13| the magnetic charge Q m is plotted against b. Delicate 
structure is observed near the critical b = b c which signifies the critical solution discussed in 
Section 5.3 (c). There may be some connection between the limiting point in the monopole 
spectrum and the critical solution. 
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Figure 10: Spectrum of black hole solutions. 

8 Stability 

It has been shown that the soliton and black hole solutions in asymptotically Minkowski and 
de Sitter space, which necessarily have at least one node in w(r), are unstable |pfl fIH| [fLOf] |j3[ . 



In contrast, the monopole and black hole solutions in the asymptotically anti-de Sitter 
space with no node in w (r) are stable for u = 0. One expects the presence of the electric 
field not to change the stability of the solitons and black hole configurations. 

In this section we give a detailed discussion for establishing the stability. We shall 
find that in asymptotically anti-de Sitter space the boundary condition for the resultant 
Schrodinger problem becomes subtle, and that the previous argument given in asymptoti- 
cally Minkowski space needs elaboration. 

8.1 Perturbation equations 

We consider small time-dependent perturbations to the static solutions to the coupled 
EYM equations. In the static solutions u{r) = w{r) = 0. In the general ansatz, ( |iT~l) and 



([L4]) we set 

u(r,t) = u(r) + 5u(r,t) 
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Figure 11: The appearance and collapse of soliton solutions in A — > 0. The bottom figure 
is the blow-up of a part of the top figure near (Qm,M) = (0, 1), in which the first four 
Bartnik-McKinnon solutions are also marked. It is seen that as A — > 0, solutions on a 
branch collapse to a point. New branches emerge as |A| becomes smaller, too. 

w(r,t) = w(r) + 5w(r,t) 

w(r,t) = 5w(r,t) 

z/(r, t) = 5u{r,t) 

p(r, t) = p(r) + 6p(r, t) 

H(r,t) = H(r)+5H(r,t) (47) 

and m{r,t) = m{r) + 5m(r,t). Substituting (pE7|) into the Yang-Mills equation (fE)| ) and 
retaining only terms linear in perturbations, one finds 

r 2 u'5p + r 2 p(8u' - 5u)] ' - (*R-*E)- ^( w 5u + 2u6w + 5tv) = (48) 

J H \ v H J H 
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Figure 12: The monopole spectrum, b v.s. M (mass) is plotted with various values of A. 
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Figure 13: The monopole spectrum, b v.s. Q M (magnetic charge) is plotted with various 
values of A. 



r 2 p(5u — 5i)) + r 2 u'5p 
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+ —(u5w + 2w5u + 28w) — — 5w 
H H 



wu 



y/5p SH\ w(l-w 2 ) r l-3w 2 r 
-i — ' 9 9 - Sp + ^ 5w = 



H \ p H J 



— (5w' + wSu) 
p 



p 



+ ^-{u 2 5w — 2u5w — w5u — 8w) 
H 



(49) 



(50) 
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uwp/Sp 5H\ (1 - w 2 ) . _ H 

—f ( — j- ) + - — = — -Sw + —w'Sv = 

H \ p H J r z p p 



The Einstein equations ( |T8"D and (|2D|) yield 



2u 2 w 2 



8m' = r 2 p 2 u\5u' — 5u) + [r 2 u 12 H — — ^JpSp + (iu 



2 2 2 
/2 P 1 U W 



H 2 



SH 



nTT , 2uwp 2 . . . •. 2w(l - w 2 ) . 
+2Hw ow H (wom + uow + ow) 

'opV 4t> f . . p 2 u 2 w 2 ( Sp SH\ p 2 uw • I 



5/7 



4vF 



— u(w'Sw — wSw') + uw 2 Su^ . 



(51) 



(52) 

(53) 
(54) 



There is residual gauge invariance specified by a gauge function fl(r, i) in (|9|) and (12). 
Making use of this freedom, one can always set either Su(r, t) — or <5z/(r, £) = 0. 



8.2 Stability analysis 

In examining time-dependent fluctuations around monopole solutions for which u(r) = 0, 



it is convenient to work in the 5u(t,r) = gauge. Eqs. (48), Q49D , and (51) become 

2pw 



(r pSv)' 



H 



Sw 



2H 



r 2 pSi) H (wSw' — w'Sw + w 2 Sp) 

p 







H 

V 



(Sw' + wSu) 



-770W + — — 

H r z p 



— Sw H w ov = (J 



whereas Eqs. (|5U|), (^), (g), and |>|) become 



#f ,/SH Sp 
— < w I 



p [ 



p 



+ <W 



p . .. w(l — w 2 ) . 

Jj sw — — S P 



1 - 3w 2 



ow = 



W = Ju/^if + 2Hw'5w' - 2wil - W ^ Sw 



Sp\ f 
P 



4^ , x , 
— u> 
r 



4v 



SH = Hw'Sw 



(55) 
(56) 

(57) 

(58) 
(59) 
(60) 
(61) 



Notice that Eqs. (53) - fl57p involve only Sv and Sw, defining the odd parity group, whereas 
Eqs. (|58|) - (|6lD involve only Sw, SH, and Sp, defining the even parity group. The number 
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of the equations is larger than the number of the unknown functions. Indeed, one equation 
in each group follows from the others. 

To derive the equation for each unknown function in a closed form, we introduce the 
tortoise radial coordinate p by 

dr = H (62) 
with which the equations for w, p, and m become 

W = ~^ W{1 ~ W) (63) 



dp 2vp / dw\ 

dp rH \dp J 



2 / Artn \ 2 

(64) 



S - +(S) ,+ ^< 1 -^} ■ (65) 

The range of p is finite, < p < p max , since p — > p and H — > |A|r 2 /3 as r — > oo: 

{r for r ~ 

3po , (66) 
Pmax - — — for T ~ OO 
I Air- 
where p = po + 0(l/r). 

In the odd parity group Eq. (|55|) expresses 5w in terms of 8v. Substituting it into Eq. 
(|56|) and making use of (|63D, one finds 

r 2 p 2 w 2 \ dp J 

5v=^P , Sw = -±±( W 0) . (67) 
r 2 p 2w dp 

Here we have supposed fluctuations to be harmonic: 5u(r,t) = e~ t<JJt 5u(p) and 5w(r,t) = 
e -*>*Sw(p). Eq. (0) follows from (||), (H), and fl63|). 

In the even parity group, (|BU|) and ( |6"T| ) express 5p and 5i7 (<5m) in terms of <5w. Eq. 
fl59|) automatically follows from flSOD , fl6"T|), (|63D and (|64|) . Eq. (|58|) becomes, with the use 
of (H), (HD and (H|, 

r rf 2 ^ 



<ip : 



2 + C^u(p) ffof = U 2 Sw 
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u, 



dp 



rH \dp J 



Again harmonic fluctuations 5w(r,t) = e 5w(p) are supposed 



(6* 



Eqs. (p7|) and (|6q) have the same form as the Schrodinger equation on a one-dimensional 
interval. Both of the potentials Up and U w are singular at p = 0, behaving as +2/p 2 . f/g 
has an additional singularity if w has a zero at p^\ Up ~ +2/(p — p^) 2 - 

The integrated energy- momentum density T ah ^J—g d 3 x due to fluctuations must remain 
finite. At the origin r = it implies that 5w = 5w = whereas 8v = 0(1). Taking advan- 
tage of the general coordinate invariance, one can impose 5p = and 5H = —2Sm/r = 0. 
At r ~ oo, Sw' , Sw' , 5v = 0(r~ 2 ). These are mild boundary conditions. One can impose 
more strict conditions such as the regularity at r = and vanishing at r = oo. As physical 
perturbations we demand that all 6w, 5w, 5u, 6H, and 6p vanish at r = oo. 

In Eq. (|67|) the potential Up(p) is positive definite. However, this does not necessarily 
mean that the eigenvalue u 2 is positive definite. It depends on the boundary condition. 



Clearly /3(0) = 0. At p 
that 



Pmax, 5v = 5w = so that (3' + h(3 — where h = w'/w. Note 



h ^ 

w dp 



|A| r 2 dw 
3p w dr 



3p w 



(69) 



where w/s are defined in 



]). For the monopole configurations with no nodes in w, h < 
{h > 0) when w is monotonically decreasing (increasing). 



Following Courant and Hilbert [28], we define 



dp {ip'(p) 2 + Up(p) ip(p) 2 } + h <p( Pn 



pmax 



dpip(p)' 



(70) 



If w(r) is nodeless, then Up(p) is regular on the interval except at p = 0. The equation 
implies that j3 = 0(p 2 ) near the origin. In this case, for an eigenfunction (3(p) in ( |B7D 
satisfies 

c 2 - (71) 

~ Jnjs) (71) 

It follows immediately that all eigenvalues u 2 are positive definite if h > so that the 
solution is stable against small odd-parity perturbations. 
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For h < more careful analysis is necessary. The lowest eigenvalue u 2 = Ai in the 
eigenvalue equation (|67D is exactly the lower bound of the set of values assumed by the 
functional V(ip, h), where <p is any function continuous on the interval [0,p max ] with piece- 
wise continuous derivatives satisfying <p(0) = 1 and M(<~p) = 1. 

\ 1 {h) = wm v V{<p]h) . (72) 

If Ai > 0, then the solution is stable against odd-parity perturbations. Suppose that <pi(p) 
saturates the lower bound for hi'. Xi(hi) = V(ipi, hi). As 

Al(/ll) = V(ipi, h 2 ) + (hi - h 2 )(fl(p max ) 2 

> A 1 (/i 2 ) + (/i 1 -/ i2 )^ 1 (p max ) 2 , (73) 

Xi(h) is a monotonically increasing function of h. Hence, if Xi(hi) > 0, then Xi(h) > for 
h > h v 

To establish the stability we utilize the residual gauge invariance. There is a zero-mode 
(with uj 2 = 0) for Eq. floTf ) with an appropriate boundary condition ho. In the A = case 
the existence of the zero-mode was utilized to prove the instability of the BK and black 
hole solutions which has at least one node in w(r) p^|||27|| . Consider the time- independent 
gauge function fl(r) in For \ fl\ <C 1, 8v — dQ/dr and 5w = —wQ. Eq. (^5|) is satisfied 
if 

^/ r _Ydn\ = 2w2Q 

dp\H dp J 



(74) 



As fi(0) = 0, Q ~ ap + 0(p 3 ) for p ~ 0. Hence Eq. ([74]) determines Q(r) up to an 
over-all constant. (3q = (r 2 p 2 /wH)(dQ/dp) is the zero mode of Eq. (|67|) and Up = Pq/Po, 
where Q' = dVt/dp etc. In this case f2(p max ) ^ and h Q = — P' /Po\ Pia ^ differs from h in 
the eigenvalue problem under consideration. If h > h Q , then Xi(h) > 0, establishing the 
stability. As d(wP )/dp = 2w 2 fl, 



, , 2w 2 Hn 

h = ho + ^r • (75) 



Nodeless solutions (w > 0) are stable if Q/Q' > at p max - 

Solving (|74|) numerically, we have determined fi(p) to find that indeed h > h for 
nodeless solutions. This analysis also shows that h becomes exactly h for the configuration 
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with w(r = oo) = 0. In this limiting case the zero mode is not normalizable; it diverges as 

(P Pmax) 

This is a general behavior. When w has a node at pk < p ma x, there appears a negative 
uj 2 mode which behaves as (p — Pk)^ 1 near p*.. 

If w(r) has n nodes, i.e. w(rj) — (j = 1, • • • , n), the potential Up develops (p — Pj)~ 2 
singularities. Volkov et al. have shown for the BK solutions in the A = case that there 
appear exactly n negative eigenmodes {uj 2 < 0) if w has n nodes |[26| . A similar conclusion 



has been obtained for black hole solutions as well p7| . Their argument needs elaboration 
in the A < case, however. 

To investigate the eigenvalue spectrum of ([J7|) in this case, it is convenient to consider 
the dual equation as was done in [2"6[|. One can write the Schrodinger equation in (|67j ) as 



Q+Q-^n = A n <p n 

Q± = ±T + f (76) 

dp (j 

where /3o(p) is the zero mode described above. <pn = 0(p 2 ) near p = and + /i<Pn = 
at Pmax- The dual equation is given by 

f d 2 ~ \ _ 
Q-Q+Vn = | - ^ + U/3jip n = \ n if n 

H 2 Sww'HQ „{w 2 HVL\'< 



^ = -^ (1 + w) -^v?^ + 8 l^wJ ■ (77) 

<p„ and <p„ are related to each other by 

f for A n + 

V" = I -1 n-1 c , n ( 78 ) 

[ <p n = /3 for A n = 
However, the boundary condition for <p n (p) depends on A n : 

<Pn(0) = 

^n(Pmax) + h n 0n(Pmax) = 

/in^-Zio-r^V ■ (79) 
n — ho 

The advantage of considering the dual equation is that the dual potential Ug(p) is 
regular except at p = where it behaves as +6/p 2 . However, the eigenvalue \ n has to be 
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Figure 14: Eigenvalues A„ (n = 1,2,3) in the dual Schrodinger equation (^) with the 
boundary condition ([79|) are displayed for monopole configurations with varying value of b. 
A = —0.01 and v — 1. When u>(oo) develops a node, the negative eigenvalue mode arises, 
signaling the instability of the solution. w(oo) = at b = 0.00168, where the numerical 
evaluation of Ai becomes difficult. 

determined self-consistently such that the boundary condition ([fjj) is satisfied. We have 
determined A n 's numerically for the monopole configurations in the lower branch in fig. 4. 
The first, second and third eigenvalues are displayed in fig. 14. One sees that Ai > for 
the nodeless configurations, but the unstable mode develops when w has a node. 

The wave function of the unstable mode in the original equation, not in the dual 
equation, diverges at the zeroes of w(r). In other words, the instability sets in around the 
zeroes of w(r). The potential Up(p) and Up(p) for the solution at b = 0.0025 are plotted in 
fig. 15. At the node of w, Up diverges, but Up remains finite. The wave function <fi(p) of 
the lowest eigenvalue (Ai = —0.0033) and the corresponding ipi(p) also have been plotted 
in fig. 15. 

For even-parity perturbation the potential U w (p) in Eq. d68|) is not positive definite. The 
first term in U w becomes negative for w 2 < 1/3. The second term also can become negative 
when w' vanishes at finite r. We have solved the Schrodinger equation (|68|) numerically 
for typical monopole solutions, and found that for the solutions with no node in w(r), the 
eigenvalues u 2 are always positive even if w(r = oo) < l/y/3. 

Hence we have established the stability of the monopole solutions with no node in w(r). 
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Figure 15: Up(p) and Up(p) for the monopole solution at b = 0.0025, A = —0.01, and 
v — 1. The wave function of the unstable mode (A = —0.0033) are also plotted. <fi(p) in 
the dual equation ([77]) is a regular smooth function, but (fx (p) in the original equation (ff6|) 
diverges at the zero of w. 

9 Summary 

New monopole, dyon, and black hole solutions to the Einstein- Yang-Mills equations have 
been found in asymptotically anti-de Sitter space. The solutions with no node in the non- 
Abelian field strengths are shown to be stable against spherically symmetric perturbations. 
The non-trivial boundary condition plays a crucial role in developing the instability for 
solutions with nodes. The stability of nodeless dyon solutions need to be established. 

Though electric and magnetic charges of monopole and dyon solutions are not quantized 
in classical theory, they are expected to be quantized in quantum theory. If this is the case, 
then at least solutions with the smallest charge would become absolutely stable. 

We have also found the critical spacetime solutions which end at finite r. These solutions 
may have connections to black hole solutions, though more detailed study is necessary. 

The solutions found in the present paper may have profound consequences in the evo- 
lution of the early universe which may have gone through the anti-de Sitter phase. We 
hope to report on these subjects in future publications. 
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